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The hydro dynamic model including the spin degree of freedom and the electro- 
magnetic field was discussed. In this derivation, we applied electromagnetism for 
macroscopic medium proposed by Minkowski. For the equation of motion of spin, 
we assumed that the hydrodynamic representation of the Pauli equation is repro- 
duced when the many-body effect is neglected. Then the spin-magnetic interaction in 
the Pauli equation was converted to a part of the magnetization. The fluid and spin 
stress tensors induced by the many-body effect were obtained by employing the alge- 
braic positivity of the entropy production in the framework of the linear irreversible 
thermodynamics, including the mixing effect of the irreversible currents. We further 
constructed the constitutive equation of the polarization and the magnetization. Our 
polarization equation is more reasonable compared to another result obtained using 
electromagnetism for macroscopic medium proposed by de Groot-Mazur. 

I. INTRODUCTION 

Hydrodynamics has been used as an effective model for discussing many-body effects and 
collective motions not only in classical systems but also in quantum ones. In the physics 
of relativistic heavy-ion collisions, it has been confirmed experimentally that the classical 
relativistic hydrodynamic model can explain the experimental behaviors of collective motions 
of microscopic quantum particles qualitatively m this case, the difference of quantum 
statistics for constituent particles of the fluid, such as boson and fermion, does not modify 
the structure of hydrodynamics. That is, it is considered that the electron fluid and the pion 
fluid obey the essentially same hydrodynamic model except for the equation of states. This 
idea is, however, not trivial under the existence of the magnetic field because of the spin- 
magnetic interaction for fermions. Then the spin degree of freedom can play an important 
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role even in constructing a hydro dynamic model coupling with the electromagnetic field. 

In fact, the spin effect under the electromagnetic field has been attracted various at- 
tentions in the physics of relativistic heavy-ion collisions 2]. There, the possibility of the 
formation of the coherent magnetic field is discussed and it is expected that the effect of 
this magnetic field to the spin degree of freedom can be observed through experiments. To 
investigate this possibility more dynamically, we need to construct a relativistic hydrody- 
namic model including the spin degree of freedom. However, such a model has not yet been 
established even in the non-relativistic case. 

For example, the spin effect is important also in the physics of quantum plasma and a 
magnetohydrodynamics with the spin degree of freedom was proposed recently |3|. In this 
model, the Pauli equation is re-expressed using classical variables such as the particle prob- 
ability density, the particle velocity and the spin vector. By assuming that each particle of a 
quantum plasma follows the same one-particle Pauli equation, the hydrodynamic equation 
is obtained by the independent sum of these constituent particles of the plasma, neglecting 
the effect of the off-diagonal quantum correlations. It is worth mentioning that the same 
re-expression of the Pauli equation with the classical hydrodynamic variables was already 



done by several groups in 1950s 



However, it is not possible to apply this result directly to the physics of the relativistic 
heavy-ion collisions, even if we ignore the difference between relativistic and non-relativistic 
models. One of the reasons is that the barotropic fluid is considered and hence the equation 
for the energy density is not derived in Ref. [3]. Furthermore, the positivity of the entropy 
production is not used in the derivation and hence the viscous tensors and the heat current 
are not obtained. Thus we need to apply the linear irreversible thermodynamics (LIT) for 
constructing a hydrodynamic model with spin. 

Another problem is that the definition of the electromagnetic momentum density in 
macroscopic medium has not yet been established js S]. In fact, there are experimental 
attempts to determine these definitions but no conclusive result is obtained because of the 
difficulty of the division of the electromagnetic part and the material part js]. 

Theoretically, there are at least three different proposals: the definitions of Minkowski 
{9], Hertz- Abraham 10] and de Groot-Mazur 11]. In the Minkowski theory, the electro- 
magnetic momentum density is defined by D x B (for the definitions of variables, see Sec. 
III). However the form of the Maxwell stress tensor is asymmetric in this theory (see Eq. 
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( 114dl) ). As is mentioned in Sec. IV B, the symmetry of this tensor is related to the angular 
momentum conservation. To introduce the symmetric Maxwell stress tensor, Abraham em- 
ployed EoHqE x H as the definition of the electromagnetic momentum density, following the 
lead of Hertz, de Groot, Mazur and Suttorp discussed this problem in two different ways, 
LIT and a certain microscopic model, and finally reproduced Hertz-Abraham's momentum 
density. However, their Maxwell stress tensor has the antisymmetric form. See Refs. jf], Q 
for details. The derived hydrodynamic models depend on the choice of these definitions. 
In Ref. js], Eu and Oppenheim derived the three hydrodynamic models using these three 
different definitions and concluded that the definition of Minkowski is promising. However, 
as was claimed in Ref. jjj, their discussion is incomplete because they could not show the 
positivity of the entropy production. 

On the other hand, Felderhof and Kroh employed the definition of de Groot-Mazur and 
succeeded in deriving the model of electromagnetic hydrodynamics which is consistent with 
the positivity of the entropy production at the last step of the derivation [3] . They introduced 
the internal angular momentum as another hydrodynamic variable, which may be interpreted 
as spin. If we cannot construct spin-hydrodynamic models consistent with LIT by using 
other definitions, the definition of de Groot-Mazur will be preferred, as was concluded in 
Ref. [7j . To confirm this conclusion, it is worth investigating once again the derivation with 
other definitions. 

In this paper, we derive the spin-electromagnetic hydrodynamics (sEMHD) by combining 
the above two approaches. First, we define the equation of motion of the spin vector so as 
to be consistent with the hydrodynamic representation of the Pauli equation. Then the 
spin-magnetic interaction in the Pauli equation induces the magnetization. On the other 
hand, the many-body effects of fluids are taken into account through the modification of the 
fluid and spin stress tensors which are obtained by employing the algebraic positivity of the 
entropy production in the framework of LIT. 

As the definitions of the momentum density, we apply the definition of Minkowski. The 
electromagnetic hydrodynamics with the definition of Minkowski was already discussed in 
Ref. fl2 |. We incorporate the information from the Pauli equation to this result. We further 
derive the explicit forms of the irreversible currents taking into account the Curie principle 
and the relaxation equations of the polarization and the magnetization, which were not 
discussed in Ref. [12 1. 
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For the purpose of the physics of relativistic heavy-ion collisions, we should discuss a 
relativistic model. However, as was discussed above, there are various unsolved problems 
even in the non-relativistic case. Thus we focus on the derivation of the non-relativistic 
sEMHD in the present work. The consistent relativistic sEMHD should reproduce the result 
of this paper in the non-relativistic limit. 

This paper is organized as follows. In Sec. II, we summarize the hydrodynamic expression 
of the Pauli equation following Refs. j^,!^]. In Sec. Ill, the Maxwell equation for macroscopic 
medium is introduced following the definition of Minkowski. Then we consider the magnetic 
polarization induced through the spin-magnetic interaction in the Pauli equation. In Sec. 
IV, we extend the result of Sec. II to hydrodynamics, so as to conserve the momentum and 
the angular momentum. To determine the stress tensors of fluid and spin, we employ the 
algebraic positivity of the entropy production, and the equation of sEMHD is obtained in 
Sec. V. The result depends on the definition of the thermodynamic pressure. By using this 
fact, it is shown that our result is consistent with that by Korteweg-Helmholtz in Sec. VI. 
The relaxation equations of the polarization and the magnetization are derived in Sec. VII. 
The comparison with other hydrodynamic models is discussed in Sec. VIII. Section IX is 
devoted to concluding remarks. 



II. HYDRODYNAMIC EXPRESSION OF PAULI EQUATION 

Let us consider the Pauli equation coupling with external gauge potentials (f> and A, 
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2m \ i c I 2 2mc 



(1) 



where c, e, m, a and -y g are the speed of light, the electric charge, the mass of a particle, 
the Pauli matrix and the g factor, respectively. The magnetic flux density B is expressed 
by using the vector potential A as B = rotA. The wave function if) has two components, 



if; 



1p2 



(2) 



Bohm et al. 0| extended the idea by Madelung [l^, Q and obtained the hydrodynamic 
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form of the Pauli equation. For this purpose, let us introduce the following three variables, 



P = V'V, 



2mi p 



A, 

mc 



(3a) 
(3b) 

(3c) 



The first two variables denote the particle probability density and the probability velocity, 
respectively. The last vector s represents the spin degree of freedom. One can easily confirm 
that the magnitude of this spin vector is normalized as s 2 = U 2 /4. 

Then all the quantities calculated from the Pauli equation can be reproduced by the 
combinations of the above three c-number variables j^, 5]. In fact, the evolution equations 
of the three variables are expressed in the closed form. Directly from the Pauli equation, 
these equations are given by 

mc 
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where d,i = di — d j. There are still terms which are expressed with the wave function ip on 
the right hand sides of the equations. To eliminate this dependence, we use the following 
relations, 



— y m 2 / e \ ' / c \k 4 v — -\ 
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(5a) 
(5b) 
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To derive these expressions, we used 
4 



2 \ n V mc 



(6a) 



-^p(sidiSj - sjdiSi) = ^2ie ijk {(ip*a k diijj) - {^a k di^*) - imps k (v l + } • 



Then the hydrodynamic expression of the Pauli equation is finally given by 
D 
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(7b) 
(7c) 



(8) 



where we introduced the material (Lagrange) derivative as 

^ = a 1 + v.v. 

The first equation represents the conservation of the particle probability. The second one 
is the equation of the probability velocity and is sometimes considered as the Euler equation 
by interpreting the first term on the right hand side as the pressure gradient % 15]- The 
last equation gives the evolution of the spin vector. That is, the above equations can be 
regarded as the ideal fluid model with spin and the external electromagnetic field. 

We will use this result to construct the evolution equation of spin in the hydrodynamic 
model. The above result is exact as the quantum dynamics of one particle, but we assume 
that the same equation is approximately applicable even to the non-interacting many-body 
systems. Then the particle probability density p is interpreted as the particle number density 
and the probability velocity v as the fluid velocity. 



III. MAXWELL EQUATION FOR MACROSCOPIC MEDIUM WITH 
SPIN-MAGNETIC INTERACTION 



Our model will be constructed so as to satisfy the total energy, momentum and angular 
momentum conservations [12J. Here total means the sum of the contributions from the 
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fluid and the electromagnetic field. For this, we have to define the energy and momentum 
densities of the electromagnetic field for macroscopic medium. In the present work, we adapt 
the Minkowski definition. 

Differently from the case of electromagnetism in vacuum, the electric and magnetic flux 
densities, D and B, are not parallel to the electric and magnetic fields, E and H, and given 
by 

D = e E + P, (9a) 
B = /i H + M + M Pmii , (9b) 

where Eq, /iq and P are the dielectric constant, the magnetic permeability and the polariza- 
tion, respectively. Here we consider two magnetizations. The quantity Mp an /j denotes the 
magnetization induced by the spin-magnetic interaction in the Pauli equation. This con- 
tribution is obtained so as to satisfy the momentum and angular momentum conservations 
as we will see soon later. On the other hand, M represents the possible higher-order spin 
contributions neglected in the Pauli equation [lsj. If we assume that the magnetization is 
induced only through the spin-magnetic interaction in the Pauli equation, M disappears. In 
this paper, this term is obtained by employing the positivity of the entropy production. See 

Sec. rvm 

As the Maxwell equation for macroscopic medium, we adapt the following equations, 

V • D = Pe , (10a) 

VB = 0, (10b) 
<9B 

VxE+ — = 0, (10c) 
<9D 

V x H — — = J, (lOd) 

where J is the electric current density, which is decomposed into the two parts: one is the 
convection current p e v, and the other the conduction current j, 

J = p e v+j. (11) 



See also the discussion of chapter XIII in Ref. 



Ill- 



Following the definitions of Minkowski, the energy and momentum conservations are 
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expressed by 

d t U ern + div S = -W, (12) 

d t g i = d k r k - f em . (13) 

Here we introduced, respectively, the electromagnetic energy density U em , the Poynting 
vector S, the electromagnetic momentum density g and the Maxwell stress tensor T lk as 

U em = -(D-E + B-H), (14a) 

S = E x H, (14b) 

g = D x B, (14c) 

T ik = E i D k + R i B k _ Uem§ ik_ ( 14d ) 

The source term of the energy density and the Lorentz force in macroscopic medium are 
given by 

W = J-E + E-d t D + H-d t B-<9 t [/ em , (15) 

f em = Pe E' + [J x B]' + YjP kQ i Ek + B k d t H k ) - d t U em , (16) 

k 

respectively 



IV. SPIN-ELECTROMAGNETIC HYDRODYNAMICS 

To obtain the evolution equation of fluids, we have to extend the result of the Pauli 
equation to the many-body case. For this purpose, we introduce new variables: the stress 
tensor of the fluid and the stress tensor of spin C y . These stress tensors are, first, fixed 
by comparing with the results in Sec. [Til an d the deviations induced by the many-body 
effects are obtained by applying LIT in the next section. 

For the discussion of hydrodynamics, it is more convenient to use the mass density defined 
by p m = where m is the mass of the constituent particle of our fluid. For the mass 
equation, we adapt the usual conservation equation, 

D . . 

J^pm = -PmV • V, (17) 

and, correspondingly, we introduce a new variable to represent the spin vector, 

s m = — . (18) 
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Then the total spin of the fluid is given by J c? 3 xp m s m . 

Of course, in the case of many-particle systems such as hydrodynamics, s m does not 
contain the complete information of spin in a many-body quantum system and our spin- 
hydrodynamic model is an approximation which will be justified when the contributions 
from the off-diagonal components of spin (interference of different spins) are absorbed into 
the many-body effect of the spin stress tensor which will be obtained by LIT. The same 
off-diagonal part is neglected even in the derivation of the spin-magnetohydrodynamics in 
Ref. J3|. As an indirect example to justify this picture, see the discussion in Sec. 9.6 of 
Ref. p], where the exact hydrodynamic expression of the two electrons system is discussed 
and the off-diagonal parts affect only a part of the spin stress tensor. In our approach, we 
assume that this off-diagonal parts is taken into account through the modification of the 
stress tensors obtained by using LIT. 

A. Momentum conservation 

We consider an arbitrary volume V whose surface is given by S. The normal vector of 
the surface is given by n. The total momentum in this region is given by the sum of the 
momenta of the fluid and of the electromagnetic field. Then the general expression of the 
momentum conservation is expressed as 

d t [ dV{p m v l + ^) = - V f dS Pm v i v k n k + V / dSp ik n k + V / dST ik n k + [ dVp m K\ 
Jv k Js k Js k Js Jv 

(19) 

Here, for the sake of generality, we introduced an external force K l which violates the 
momentum conservation. This should be satisfied for any volume and then we obtain the 
following equation, 

p^w/ = S d ^ + pmK ' L + £»• (20) 

j 

To reproduce Eq. (1751) . we identify 

K* = 0, (21a) 
P ij = (PT + (PT, (21b) 
(p-)« = {ff + lnp m - p m J2(diS l J(d jS l J + ^^s m ■ B, (21c) 
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where (p s )^ and (p a ) iJ are the symmetric and anti-symmetric parts of the stress tensor of 
the fluid for the exchange of the indices i and j. The terms (p s )^ and (p 11 )^ represent the 
modifications of the stress tensor due to the many-body effect which will be determined later 
by employing LIT. 

The term Ylk s m^i^ k in Eq. (17b 11 is reproduced from the magnetization in the Lorentz 
force. In fact, by choosing 

(22) 

the Lorentz force is expressed as 

f em = p e E l + [JxBf + V {D'd.E 1 + B%-(B - M) 1 ) - diU em - 1iBiy^B l dd m . 

i I f* J 2c i 

(23) 

Then, by combining the contributions from the last terms of Eqs. (I21c[) and (123]) . the term 
J2k s m diB k is reproduced. 
In short, we can construct the equation consistent with the momentum conservation when 
we choose Nlp au u by Eq. ( 1221) . The consistency of this Nlp au u is confirmed in the discussion 
of the next subsection, the angular momentum conservation. 

B. Angular momentum conservation 

Usually, in particular in the relativistic case, hydrodynamic models are derived so as to be 
consistent with the conservation of energy and momentum but the angular momentum is not 
considered explicitly, because this conservation is automatically satisfied when the energy- 
momentum tensor is chosen to be symmetric. In the present case, however, the Maxwell 
stress tensor is not symmetric T*- 7 ^ and we need to employ the angular momentum 
conservation for the construction of the consistent model. 

The total angular momentum density is given by the sum of the fluid part r x p m v, the 
electromagnetic part r x g and the spin part p m s m . Then the conservation for the angular 
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momentum of the total system is expressed as 



d t / dV[r x (p m v + g)f + d t dVpmsi 
Jv Jv 




where G, if it exists, means the source of the angular momentum. That is, the last term on 
the right hand side represents the violation of the angular momentum conservation. From 
this expression, we obtain the following differential equation, 

P™§-/m = E € ^P kj + E ^ + P x E + B x H] 1 + p m G\ (25) 

jk 3 

To derive this expression, we used the momentum equation obtained in the previous sub- 
section. Note that the third term on the right hand side of Eq. ( 1251) disappears in vacuum 
because D || E and B || H. 

Using the definition of Mp mK given by Eq. (1221 . Eq. (!7c|) is reproduced by the following 
identification, 

G i = 0, (26a) 
C ij = Pm [s m xd j s m \ i + &\ (26b) 

where the possible modification of the spin stress tensor by the many-body effect is denoted 
by C y . The interaction term s m x B in Eq. (ITcl) is induced from B x H in Eq. ( |25l) through 

^-Pauli- 

V. LOCAL THERMAL EQUILIBRIUM AND ENERGY EQUATION 

In this section, we calculate the modifications of the stress tensors, (p s ) lJ and C u by 
employing LIT. 

For this purpose, we first calculate the energy equation. Suppose that spin can be ex- 
pressed as s m = Id using the moment of inertia / per unit mass and the angular velocity O 
as in the case of the classical angular momentum. Then the total energy density is expressed 
as 

P^(v 2 + m 2 )+p m U + U em . (27) 
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These four terms represent, respectively, the kinetic, rotational, internal and electromagnetic 
energy densities. Here we introduced the internal energy per unit mass by U. The validity 
of the assumption of the spin rotational energy used here should be argued carefully. We 
will discuss this point in the concluding remarks. 
The energy conservation is, then, expressed as 



d t / dV 



^(v 2 + /0 2 ) + p m U + U e , 



v ■ n) + / dS ^ P 



lJ n J v l 



dS \^ + m 2 )+ Pm u 

+ / dSJ2c ij n j Q i + [ dV Pm {K-v)+ [ dV Pm (G-n) 
Js {j Jv Jv 

- [ dS(S-n) - [ dV(V-q), 
Js Jv 

where q denotes the heat current. The corresponding differential equation is given by 

j • [E + (v x B)] - 



(28) 



Pm Dt U Dt Uem 



D* D* 
D — E + B — H 

Dt Dt 



(29) 



ijk ik ik 

where we introduced the material derivative with rotation as 

§lx = ^x-nxx, (so) 

for an arbitrary vector X. 

In general, (p s )* J has 6 degrees of freedom. In this work, following the usual phenomeno- 
logical derivation of the Navier-Stokes equation, we parametrize it by the three quantities, 
the thermodynamic pressure P, the bulk viscous pressure t d and the shear stress tensor 
(t s )^', 

h 2 



(PT + T-zPmdA lnp m - p m J2(djS l J(d k s l J = -P5» + r D 5 l3 + {r s f 



4m 2 



(31) 



where (r s ) ij = (r s ) ji and EiO 5 )" = °- The term S?Pm,djd k In p m - p m EK 9 i s m)( 9 fc s L) is 
known as a kind of pressure comes from quantum fluctuation 4j, [5] . In the above expression, 
we define the pressures and the stress tensor so as to renormalize this quantum effect into 
their definitions. It is possible to consider the more complex decomposition of the fluid 
stress tensor. See p. 311 of Ref. [ll| for details. 



13 



Substituting p*- 7 and C iJ into Eq. (129|) . we obtain the energy equation, 

D TT D „ D* D* „ 1 P> 

Dt Dt Dt Dt p m Dt 



+j ■ [E + (v x B)] + [r D + ^(s m ■ B)j 9 + $>^ fc ^' fc 

-p2 • (n - + ^ c jk d k n j - v • q , (32) 



where 



y-> i)jr' - i),r J . (33b) 
(^f^-^e.y", (33c) 

jk 

= djV 1 + dy - (33d) 
e = diV\ (33e) 

i 

In order to obtain the thermodynamic relation of our system, we introduce the following 
quantities per unit mass, 

Uem = (34a) 

Pm 

D = — (34b) 

Pm 

B = — . (34c) 

Pm 

By rewriting the energy equation with these variables, we find that the total energy, which is 
given by the sum of the internal energy per unit mass U and the electromagnetic energy per 
unit mass U em , can be regarded as a function of p m , D, B and additional contributions which 
will be interpreted as heat. From this behavior, we assume the following thermodynamic 
relation in this work, 

dU + dU em = E • dt) + H • dB - (P + U em )d— + TdS en , (35) 

Pm 

where T and S en are temperature and the entropy per unit mass, respectively. The similar 
thermodynamic relation is obtained in Sec. 5 of Ref. 



. This result means that the 



thermodynamic work is modified from Pd(l/p m ) to (P + U em )d(l/ p m ) — E-dD — H-c£B. This 
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can be shown by calculating the thermodynamic work in the quasi-static process explicitly, 
as is discussed in Appendix |X] Note that the contributions from the electromagnetic field 
appearing in the above relation cancel and the usual thermodynamic relation is reproduced 
when the polarization and the magnetization disappear, P = M + ~M.p au u = 0. 
Then the evolution equation of the entropy per unit mass is given by 



,DS K 



DU DU P 



+ D 



D*E 



D*H 



Dl/Pr 



, B ■ ——— + (P + f era )— (36) 

Dt Dt Dt Dt Dt K em) Dt 1 ; 

Here we interpret D/Dt as D* / Dt when it is operated to vectors. In fact, there is the 
following relation, 



D D* D* 

X Y = X — Y + X • — Y, 

1-S b 



Dt Dt 
for arbitrary vectors X and Y. 

This equation is cast into the following form, 



(37) 



D 



S P 



Pm "en 

where the entropy production is defined by 



-V- 



q 

T 



(38) 



0\s 



ij.(E + vxB) 



jk 

■i P ; ■ (o - IcA + ^c ik d k w + q • 

^ ' ik 



(39) 



The unknown currents and tensors are determined by employing the algebraic positivity of 



this entropy production 



11] . Then we obtain 



j = cr(E + v x B 

D 



T 



(6 - ■ P, 



2c 



.S\ik 



Jk 



t ) = r)e 
D <9 V 



q = 



l 

«v-, 



-a n 



i 



(40a) 
(40b) 
(40c) 
(40d) 
(40e) 

(40f) 



The first equation corresponds to the Ohm law. The electric conductivity, the spin diffusion 
coefficient, the bulk viscosity, the shear viscosity and the heat conductivity are given by a, 
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D s , £, rj and k, respectively. The last equation represents the contribution from spin and 
the vorticity of the fluid, and its transport coefficient is given by A. 

From LIT and the Curie (symmetry) principle 11[ , the most general irreversible flows are 
given by the linear combinations of currents which have the same transformation properties. 
For example, both of J and q are vector currents and the most general expressions of these 
two currents are given by the linear combinations, 



r 1 E + v x B 

<lt0t — -^11 V— + — ; 

1 E + v x B 

hot — -t>2lV— +L 2 2 — • 



(41a) 
(41b) 



The coupling between the charge and heat currents are known as the Seebeck and Peltier 
effects. The constraint to L12 and L21 is given by Onsager's reciprocal theorem 111 ]. 

We did not consider the separation of the symmetric and anti-symmetric parts of C'^ as 
was done for (p s ) u . If we consider such a separation, it is, in principle, possible to consider 
the mixtures between djS l m and djV 1 as was done for q and j. This mixture is, however, not 
considered here. See also the discussion in Sec. HXl 

Finally, the equations of the spin-electromagnetic hydrodynamics are summarized as 



D 

Di Pm 
D i 

Pm V 

h Dt 



"Pm(V • V), 

-d i (p-ce) + Y / d j ( V e 



(42a) 



j 



-p e [E + v x Bf + [ Mot xBf + J2 [ Dl 9iE l + B'd.H 1 ] - d^ e . 



(42b) 



D 
Dt 



D 

Pm Di U 



T 



-UJr. 



x djS m }i + D s djS l m ] 



— D x E + — B x (M + M Pauli 
£ o Po 



(42c) 



§-U em + E ■ + H ■ ^B + j« • [E + (v x B) ] + P±-^ Pm 

+ \ 5>' fc ) 2 + A (t " 2 + t 5><) 2 - v ■ ( 42d ) 
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VI. EQUATION OF STATE DEPENDENCE 



In this section, we consider a special case where the the properties of the polarization 
and the magnetization are restricted as follows, 



P 

M + M Pmli 



(43a) 
(43b) 



where xe and xh are the electric and magnetic susceptibilities, respectively. That is, D (B) 
is proportional to E (H). Then the thermodynamic relation (135]) is reduced to 



TdS P 



dU + {P + U em )d—-\ 

Pm £ 



E 2 d^ + H 2 d^ 



(44) 

Pm Pm _ 

So far, we have used the thermodynamic relation among dU, dS en , dp m , d(xE/p m ) and 
d{xn/pm)- For the later convenience, we choose xe and xh to be thermodynamic variables 
instead of Xs/Pm and Xn/Pm, respectively. Then the thermodynamic relation (142d|) is, 
again, re-expressed as 



TdS e 

nn 



dU + Pd— - [E 2 d XE + H 2 dx H ] 

Pm ^Pm 



(45) 



Following Refs. [111 . Il2l Il8| , we consider the free energy F = U — TS en and its variation, 



dF = d(U - TS en 
P + p 



2 dp r , 



+ Pr 



H 2 d X H 

2 dp m 



Pm 



<? — 



E 2 d XE H 2 d X H 



2p m &T 2p m &T 



dT. 



(46) 



Here we assumed that xe and xh are functions of p m and T. Thus we can eliminate the 
Xe and xh dependences by introducing the following new pressure and entropy which are 
functions only of p m and T, 



dF = -P d— - S dT, 

Pm 



(47) 



where 



P = P (pm,T) - p., 



&dxi 
2 dp n 



H 2 d x 



H 



2 dp r , 



E 2 &xe H 2 &xh 

Sen = So (p m , T) + - — — - + - — — . 

2p m dT 2p m dT 



a) 
(48b) 
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By using this new pressure, our momentum equation is re-expressed as 

Pm Wt vt = ^ (P ° ~ C6) + S 9 i^ e<i ) + Pe[E + v x B]« + [j tot x BY 

~T d * XE ~ ~Y daH + di [ pm Yd^ ) + di [ pm ^d^ ) ■ (49) 

Here, for the sake of comparison, we ignore the contributions from spin. The second line on 
the right hand side is equivalent to the so-called Korteweg-Helmholtz force, which is known 
as one of the possible candidates for the interactions between fluid and the electromagnetic 
field [17] . The Kelvin force is another candidate and the difference between the Kelvin and 
Korteweg-Helmholtz forces is absorbed into the definition of the pressure again. See the 

discussion in Refs. HQ. 



VII. DYNAMICS OF POLARIZATION AND MAGNETIZATION 



We have considered that the electromagnetic field is the solution of the Maxwell equation 
for macroscopic medium. On the other hand, as is discussed in Refs. [7, s[ [ll], it is also 
possible to construct constitutive equations of the polarization and the magnetization by 
introducing the two different electromagnetic fields: one is the solution of the Maxwell equa- 
tion for macroscopic medium and the other the conjugate variables of the polarization and 
the magnetization appearing in the thermodynamic relation. The latter fields are denoted 
by E e? and H eq , and then the thermodynamic relation is re-expressed as 



dU+—d-CE-P 

Pm 2 



H M 



tot) 



E eq ■ dP + W q ■ dm 



lot 



Pd h TdS en , 

Pm 



(50) 



where M tot = M + M Pauli . 

By using this relation, we can calculate the entropy production again and obtain 



-j-(E + vxB) 



+ u D + ^(s m -B; 



2c 
CI - 



ik 



4 



T 2T 

fC* k d k W + q ■ V 

ik 



1 

T 



1 D* 1 D* 

+-(E - E eq ) — P + -(H - W q ) — M 

T v ' Dt T K ' Dt 



tot 



(51) 
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The new contribution to Eq. ( )39|) comes from the last two terms on the right hand side 
of the equation. The positivity of the entropy production for these terms then leads to 

— P -fixP = — (E - E eq ), (52a) 
Dt T p 

^-M tot - ft x M tot = — (H - H e «), (52b) 
Dt t m 

where r v and r m are new transport coefficients, characterizing the time scales of the relax- 
ations. 

There are several proposals for the evolution equations of the polarization and the mag- 



netization. In Refs. |8j, |ll|], the derived equations of the polarization and the magnetization 
do not contain the terms representing the rotation of the polarization and the magnetiza- 
tion, which corresponds to the second terms on the left hand side of the above equations. 
This is because the internal angular momentum or spin is not considered as a hydrodynamic 
variable. 

The same rotation terms are obtained in Ref. [3] where the rotational vector ft is replaced 
by (V x v)/2, and it is discussed that the derived equation of the polarization is close to 



Eq. (2.6) of Ref. 19J1. If we can use the same replacement in our result, one can see that 

n ft 

our equation is more similar to Eq. (2.6) of Ref. [19] compared to that of Ref. [7|. 

VIII. COMPARISON WITH OTHER THEORIES 

As was mentioned in the introduction, there are three different hydrodynamic models 
including the spin degree of freedom (internal angular momentum) and the electromagnetic 
field. In this section, we would like to summarize the comparison among them. 

n 

In Ref. [3(, the spin-magnetohydrodynamics is obtained by using the result by Bohm et 
al. As was mentioned in the introduction, however, the equation of the energy density is 
not obtained. Moreover, the form of the stress tensor, which corresponds to viscosity, is not 
calculated in the framework of this theory. 

In Ref. the model is obtained phenomenologically using the definition of de Groot- 
Mazur. In this derivation, the positivity of the entropy production is satisfied in the last of 
the derivation using LIT. However, the form of the fluid stress tensor is assumed to be given 
by Eq. (4.7) of Ref. [7| and the approach of LIT is not completely employed. Because of 
this incompleteness, the mixing effect of irreversible currents such as the Seebeck and Peltier 
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effect is not discussed. Furthermore, the internal angular momentum included in this theory 
is not equivalent to spin. For example, the traditional conservation equation is assumed for 
the internal angular momentum, 

d t s v + V • (s„v) = source terms, (53) 

where s„ is the internal angular momentum density of Ref. Q. See Eq. (4.3) of Ref. 
0]. However, this structure is different from our spin equation (142c| . which is the non-linear 
equation of s m . Thus it is not adequate to identify the internal angular momentum discussed 
in Ref. {t| with spin. Moreover, as was shown in the previous section, our polarization 
equation is more promising than that of Ref. [?| ■ 

Our approach is the extension of the argument of Ref. [12j to the case with the spin 
degree of freedom as was mentioned in the introduction. Thus same thermodynamic relation 
and electromagnetism for macroscopic medium are used. However, the internal angular 
momentum introduced in this work is essential same as that of Ref. |7| and hence is different 
from our spin equation. Moreover, the mixing effect of the irreversible flows and the dynamics 
of the polarization and the magnetization are not investigated. 

If the spin degree of freedom (or the internal angular momentum) is neglected, sEMHD is 
reduced to electromagnetic hydrodynamics. For the purpose of reference, we will comment 
some results of this ty pe o f theory. 

The result of Ref. 



ll| is the essentially same as Ref. [7| except for the absence of the 



a 

internal angular momentum. However, differently from Ref. [7J, the fluid stress tensor is 
determined by employing the positivity of the entropy production completely as was done in 
this paper. Moreover, the mixing effect of the irreversible currents is also taken into account. 

n 

In Ref. [8J, the derivations of electromagnetic hydrodynamics are considered for three dif- 
ferent choices of electromagnetism for macroscopic medium: the Minkowski, Hertz- Abraham 
and de Groot-Mazur definitions. It is concluded that the choice of Minkowski is promising, 
but the derived entropy production does not satisfy the algebraic positivity and the form 
of the viscous tensor and the heat current are not evaluated. The thermodynamic relation 
used in Ref. [8[ (given by Eq. (4.17)) is different from what we used in this work and this 
is the reason why the positivity of the entropy production is successfully satisfied in our 
derivation. 



The thermodynamic relation used in Ref. 



18} is the same as ours (See Eq. (5.10b) in Ref. 
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18||), while the Abraham-Becker stress tensor is used for electromagnetism for macroscopic 
medium. Moreover, the form of the fluid stress tensor is assumed and is not obtained by 
using the positivity of the entropy production. 

IX. CONCLUDING REMARKS 

The hydrodynamic model including the spin degree of freedom and the electromagnetic 
field was derived. In this derivation, we applied electromagnetism for macroscopic medium 
proposed by Minkowski. For the evolution equation of spin, we assumed that the hydrody- 
namic equation of the Pauli equation is reproduced when the many-body effect is neglected. 
The many-body effect for the fluid and spin stress tensors were obtained by employing the 
algebraic positivity of the entropy production in the linear irreversible thermodynamics. 

As was shown in this paper, the positivity of the entropy production is satisfied in the 
derivation of sEMHD even if we use the definition of Minkowski. This conclusion is different 



from that of Ref. 



a, 



where it is claimed that the definition of de Groot-Mazur should be 



employed for the consistent theory. In fact, our derivation is consistent with LIT more than 
Ref. |7|, because the assumed form of the stress tensor is not employed in our approach. 



Furthermore, our polarization equation can reproduce the result of Ref. 19[ better than 
Ref. IJ. 

However, it does not necessarily mean that we should choose the Minkowski definition, 
because it may be still possible to derive a consistent hydrodynamic model in the de Groot- 
Mazur definition by considering the different thermodynamic relation. Similarly, it is also 
important to confirm the applicability of the definition by Hertz- Abraham to sEMHD. The 
superiority of these derived hydrodynamic models will be studied by, for example, the sta- 



bility analysis of the models as was done for the relativistic hydrodynamic models 20]. 

In the present argument, we discussed the energy conservation by considering the spin 
vector per unit mass as an analogue of the classical angular momentum. Then the spin 
rotation energy is assumed to be given by JO 2 . In fact, it is known that the kinetic energy 
of a quantum particle with spin can be expressed by the sum of the kinetic term of velocity 



and the rotation of spin, as is shown in Ref. 21] 
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where 

= m 

((VlnV*V) 2 )' 

Here ( ) represents the quantum mechanical expectation value, respectively. One can see that 
this I is the essentially same as the classical moment of inertia. In fact, when the Gaussian 
distribution is considered, ip*ip ~ e~ r2 ^ 2a2 \ we obtain I = ma 2 . Thus our classical treatment 
of the spin degree of freedom can be justified when the fluid is given by the ensemble of 
wave packets of constituent particles of a fluid with a well-defined width a. 

We developed the model of sEMHD by taking notice of the classical aspect of the spin 
degree of freedom. Then the many-body behaviors of the original quantum dynamics is 
not reproduced exactly. For example, the effect of the Pauli exclusion principle will not 
be included. Thus our model is applicable for high temperature systems where this effect 
becomes relatively small. 

We considered the mixture of the two vector currents and obtained the Seebeck and 
Peltier effects. In principle, it is also possible to consider the similar mixture for tensors. 
For example, the shear stress tensor can be modified by the mixture with the spin tensor as 

ije* — ► L£e« + Lge?, (56) 

with 

e? = djsl - d^. (57) 

The magnitude of the mixture is determined by Onsager's reciprocal theorem. 

In this discussion, we did not consider the multi-component fluids with different charges. 
The generalization of our approach to such a case will be implemented following the argument 



ofRef. [11 1 . 



As was discussed in the introduction, we would like to construct the relativistic hydrody- 
namic model including the effect of spin, but the generalization of the present argument to 
the relativistic case is not trivial, because spin is not an appropriate quantum number for 
relativistic systems. However, the consistent relativistic theory should reproduce our result 
in the non-relativistic limit, as the relativistic hydrodynamic model reproduces the Navier- 
Stokes equation in the appropriate non-relativistic limit. As one of examples of relativistic 



hydrodynamics associated with spin, see Ref. 22] 
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Appendix A: Thermodynamic Relation 



In this appendix, we derive the thermodynamic relation given by Eq. fl35|) . 
The general expression of the first law of thermodynamics is 

AQ = dE + dW, 



(Al) 



where AQ and dW represent heat and the thermodynamic work, respectively. These quan- 
tities are expressed by thermodynamic variables in the quasi-static process. In particular, 
the usual thermodynamic work in this process is expressed as 



dW = PdV. 



(A2) 



When we have the electromagnetic field, it acts as an external field to a thermodynamic 
system and gives rise to another contribution to the above thermodynamic work 23j. The 
contribution is estimated as follows. 

We separate a fluid into the ensemble of fluid elements (cells) and assume that the thermal 
equilibrium is satisfied inside of these elements. On the other hand, the change of the 
momentum of the fluid is induced by the stress tensors acting on the surface of the fluid 
element. Then the change of the work for a fluid element is defined by 



Dt ^ 



ik 



dSv\p ik + T lk )n k . 



ik\ 



This quantity is calculated as 
D 



Dt 



-W 



+ 



dV 
dV 



D 



(A3) 



j • (E + v x B) + p m E— D 



p m H— B 

Hm Dt 



(A4) 



To obtain the correction to the thermodynamic work PdV (or Pd(l/p m ) in the notation 
of hydrodynamics), we consider the quasi-static process which is reversible. Then we ignore 
the contribution corresponding to the higher order of the fluid velocity, v 2 = and the 
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contribution from the irreversible flows. The thermodynamic work in the quasi-static process 
is finally obtained by 



dW 



dVp r , 



-{P + U em )d— + ErfD + HdB 



Pi, 



(A5) 



Thus one can see that, because of the spin (the magnetization) and the electromagnetic 
degrees of freedom, the thermodynamic work is modified as 



Pd- ► (P + U em )d— - Edt) - HdB. 

Pm Pm 

As a result, the thermodynamic relation is modified as Eq. fl3"o"|) in the present case. 
The similar result is obtained in Ref. 18] . 
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